Dissipative solitons stabilized by a quantum Zeno-like effect 
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The solitary wave 1 (or the soliton 2 ) is a ubiquitous nonlin- 
ear phenomenon, which has been widely observed in many 
systems, particularly, in optical fibers 3,4 and Bose-Einstein 
condensates (BEC) 5,6,7 . It was recognized that the stability 
of a soliton is due to a subtle balance between competing 
features like dispersion and nonlinearity in a medium. The 
practical applications of the soliton raise the requirement 
for optimal soliton management 8,9,10 . However, due to the 
presence of dissipation, the soliton inevitably fades away. 
Here we show that the dissipative effect can be completely 
suppressed by an analogue of the quantum Zeno effect 11 (in 
short we will call it quantum Zeno-like effect) and propose a 
novel way to stabilize the soliton. The proposal opens a new 
perspective to realize an ideal optical soliton transmission, 
a seemingly impossible dream in communication technolo- 
gies. It is also of fundamental interest to the quantum Zeno 
effect itself and the matter-wave dynamics. 

In quantum mechanics, the quantum Zeno effect shows 
that an unstable particle can be completely stabilized 
by frequent measurements 11 , which has been confirmed 
by many experiments 12,13,14 and its possible applications 



have been extensively explored 15,16 . Solitons are local- 
ized, finite energy states in a medium. They behave like 
elementary particles, can pass through each other and 
preserve their shapes and speeds after collisions. How- 
ever, the presence of dissipation effect makes the soliton 
behave like an unstable particle, so it will inevitably de- 
cay and eventually disappears in a medium. The insight 
gained from the quantum Zeno effect in quantum me- 
chanics provides a possible mechanism to stabilize such 
a soliton if it is frequently "measured" . Through a rigor- 
ous mathematical analysis, the soliton is, indeed, shown 
to become stable if the external potential applied and/or 
the Feshbach resonance along with the dispersion is pe- 
riodically modulated with high frequency. This is due to 
the fact that the "inevitable" dissipation effect is com- 
pletely suppressed by such a frequent "measurement" , 
even in the presence of the dissipation in reality. 

We study the dynamics of the soliton governed by 
a one-dimensional (ID) dimcnsionless nonautonomous 
nonlinear Schrodingcr (NLS) equation 
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Here u(x, t) can represent either the envelope of pulses 
in nonlinear optics with different notations or the macro- 
scopic wave function of BEC with strong transverse con- 
finement (i.e., effective ID model) at the mean-field level. 
Therefore, the soliton satisfying equation (1) can be the 
temporal (if V(t) = 0) or the spatial optical soliton or 
the matter-wave soliton. For convenience, below we use 
the notations for the matter-wave solitons. f(x,t) and 
g(x,t) arc dimcnsionless and denote the time- and space- 
dependent dispersion and nonlinearity managements, re- 
spectively. Here the space- and time-coordinates x and 
t are measured in units of a r and lj~ , where a r and ui r 
being the transverse harmonic oscillator length and the 
transverse confining frequency of BEC, respectively, e 
and S are constants in the standard NLS equation (see 
Method). V(t) represents the external harmonic poten- 
tial applied, whose strength and sign can be tuned ex- 



perimentally. The linear dissipation (7 < 0) or gain 
(7 > 0) rate is also time-dependent. Generally, these co- 
efficients should be real. Due to the introduction of these 
control parameters, the soliton obeying equation (1) can 
be called as nonautonomous 17 or dissipative 18 soliton or 
similariton in nonlinear optics 19,20,21 . 

The integrability study of the nonautonomous N1S 
equation has a long history, see, for example, ref. (9). A 
complete integrability condition of equation (1) neglect- 
ing dissipation and/or gain (i.e., 7(f) = 0) was clearly 
shown in ref. (17) and further studied recently by the 
Painleve analysis 22,23,24 . However, it is nontrivial to ex- 
plore the integrability of equation (1) in the presence of 
the dissipation and/or gain effect. Below we show a gen- 
eralized "integrability" condition of equation (1). It can 
be obtained by the Painleve analysis 22 . For details of the 
Painleve analysis one can refer to ref. (23). It reads 
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Here we would like to point out that the dispersion and 
nonlinearity managements arc not allowed to be space- 
dependent, even if their space-dependence assumptions 
are initially made in equation (1). This is a necessary 
condition for equation (1) to pass through the Painlcvc 
test. Hereafter the coordinate x in f(x,t) and g{x,t) 
is omitted. In addition, we also mention that equation 
(2) is essentially not a complete integrability condition in 
a strict physical sense, namely, the infinite conservation 
laws are not assumed. It is the condition under which the 
Painleve analysis can pass through. In this sense, we will 
call it the Painlcvc integrability condition, which sheds a 
new light on the "integrability" behavior of equation (1). 

The seemingly complicated Painleve integrability con- 
dition obtained for the first time shows a subtle balance 
between four parameters f(t),g(t),V(t) and 7(f) affect- 
ing the dynamics of the nonautonomous soliton, and thus 
has a profound implication to the control of the soliton 
dynamics. Any three parameters out of them can be set 
independently, while the remaining one can be tuned ac- 
cording to equation (2). In this case, the soliton satisfied 
equation (1) can be managed as a whole in the sense of 
the Painlcvc integrability. This feature provides many 
possibilities to control the soliton dynamics in different 
fields such as nonlinear optics and BEC, deserving fur- 
ther exploration. 

Here we focus on the dynamical stability of the soliton, 
which is of fundamental importance to the study of the 
matter-wave dynamics and in particular to the propaga- 
tion of the optical solitons in transmission lines. In this 
case, f(t),g(t) and/or V(t) are set independently and 
7(t) is effectively determined by these three independent 
parameters. At first glance, this assumption seems to 
be logically incorrect since the dissipation and/or gain 
mechanism always physically exists in reality. However, 
when some control parameters are modulated, the ac- 
tual 7(t) acting on solitons might be different from that 
without such modulations. Here it is quite instructive to 
think of a comparison with the quantum Zeno effect. A 
physically unstable particle always decays with a natu- 
ral decay rate named as Rq. However, it is well known 
that that frequent measurements, for example, with a 
frequency u>, on the unstable particle can slow down its 
decay. This fact tells us that the actual dissipation rate 
named as of the unstable particle with frequent mea- 
surements must be different from Rq. In this case, Rq is 
only a nominal, while R u is the actual dissipation rate 
of the unstable particle. The soliton case is similar. The 
only difference is that the soliton is not a real particle, 



though it behaves like an elementary particle. For con- 
venience, we take the dissipation and/or gain as 7 Ii) (t) 
when dispersion, nonlinearity or the external potential is 
modulated periodically with frequency u. 

Focusing on the dynamical stability of the soliton, we 
rewrite equation (2) as 

(<) - ll (t) - a(t)^ (t) + 0(t) = 0, (3) 

where a(t) = jf^f t f(t) - -^jf t g{t) and (3(t) represents 
the remaining terms in (2) independent of J u (t). This 
is a standard Riccati equation 25 . Mathematically, the 
Riccati equation can not be solved analytically in a gen- 
eral case. Here we resort to the Rungc-Kutta method 
to solve numerically this equation. We first consider the 
effect of the external potential V(t) by fixing f(t) and 
g(t) to f{t) = git) = ±1 and let e = 1/2 as usual. Thus, 
a{t) = and /3(t) = TV it). 

Figure la shows 7^(i) as a function of time. At low fre- 
quencies, 7n;(i) oscillates and turns slowly up with time. 
Finally, at certain moment Juj(t) diverges. With increas- 
ing frequency, the diverging point of 7 w (t) is postponed 
rapidly and finally, at the time interval wc studied (up 
to 100 here and not fully shown for clarity) no diver- 
gence is observed at least for u> > 128. Meanwhile, the 
magnitude of oscillation decreases quickly. This clearly 
indicates that the matter-wave soliton can be stabilized 
by applying a modulated external potential with high fre- 
quency. This can also be seen from the behavior of T w (t) 
with increasing frequencies, as shown in Figure lb. When 
lu = 2048, no visible deviation from zero is observed for 
the whole time domain studied. The explicit time evo- 
lution of the nonautonomous bright soliton is shown in 
Figure le-f with frequencies of u> = 2,8, and 32. With 
increasing frequencies, the nonautonomous bright soli- 
ton approaches gradually to the canonical bright soliton 
shown in Figure Id. The overlap integral of the evolv- 
ing bright soliton with its initial state is shown in Figure 
lc for different frequencies. At high frequencies the soli- 
ton evolves slowly and the deviation from its initial state 
keeps invisible. Analytically, the nonautonomous soliton 
can be effectively reduced to the canonical soliton form 
since r u — > at the limit of high frequency. As a result, 
a(x,t) = 0,X(x,t) = x,T{t) = i, andc(i) = neglect- 
ing the constants C\,C2, and C3 (see Method). In this 
case, the Painleve integrability condition recovers con- 
sistently to the complete integrability condition 17,24 . As 
a result, the nonautonomous soliton becomes stable and 
completely integrable. This is a quantum Zeno-like effect 
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Here u{x,i) obeys equa- 



FIG. 1: The effective dissipation and/or gain rates and 
the time evolution of the nonautonomous solitons in 
the periodically modulated external potentials, a, The 

effective dissipation and/or gain rate 7w(t)- b, The effective 
accumulated dissipation and/or gain rate, which is defined 
by T^(t) = J* ~fuj(t')dt' . c, The normalized overlap integral 

J u* (x ,0)u(x ,t)dx 
J |u(a:,0)| 2 d^ 

tion (1) and u(x,0) is its initial state, which is equivalent to 
the canonical bright soliton. d, The canonical bright soliton 
Q{X, T) = sech(X) exp(iT/2). e - f, The explicit evolution of 
the nonautonomous bright soliton in periodically modulated 
harmonic external potentials with increasing frequencies from 
uj — 2 to 32. The time interval for the soliton evolution is 
from to 10 (4 in the case of u = 2). Here the frequency 
and time are measured in units of uj r and av* , respectively. 
The strength of the harmonic external potential is modulated 
by V(t) = V cos(wt) with V = 1. Other parameters used are 
f(t) — g(t) = 1, e = 1/2 and for all cases the initial value 
of "/u(t) is set to be zero. The result of f(t) — g(t) = —1 is 
similar. 
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FIG. 2: The effective dissipation or gain rates and the 
time evolution of the nonautonomous solitons with 
periodically modulated nonlinearity. a - g fully corre- 
spond to those in Figure 1. The time interval for the soli- 
ton evolution in e-g is plotted from to 10. The nonlin- 
earity modulation is given by g(t) = exp[Sa cos(urf)], where 
5a — 1/lo 2 is taken. Other physical variables used are 
V(t) = and f(t) = 1. 



in the context of nonlinear dynamics found here for the 
first time. It builds up a clear relationship between the 
nonautonomous and the canonical solitons. This surpris- 
ing result unambiguously shows the particle nature of the 
nonautonomous solitons and provides a dynamical way to 
stabilize them. 

Experimentally, it is easy to modulate the external po- 
tential by changing periodically the direction and mag- 
nitude of the magnetic field controlling the magnetic 
traps for BEC. Actually, an easier way to stabilize the 
matter-wave solitons is to use the Feshbach resonance 
management 26 . To show this, let f(t) = 1 and g(t) = 
exp(Sacos(tut)), where da and uj are the amplitude and 
the frequency of the Feshbach resonance modulation. For 
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simplicity, we also take V(t) = 0. In this case, one has 

a(t) = — 25a lu sm(ujt), 

(3{t) = -8a uj 2 (cos(wi) - 5a sin 2 (cot)) . 

Limiting to weak disturbances, we take Sa = 1/uj 2 . The 
result is shown in Figure 2a & b. Similarly, the matter- 
wave soliton can be stabilized by a high frequency mod- 
ulation Fcshbach resonance. As an example, Figure 2 e-f 
shows the nonautonomous bright solitons evolving with 
time. Likewise, with increasing frequencies, the nonau- 
tonomous soliton approaches gradually to the canonical 
soliton evolution, as shown in Figure 2d. For uj = 2048, 
the nonautonomous soliton has no visible deviation from 
its initial state, as shown in Figure 2c. This is another 
example that the nonautonomous soliton can be stabi- 
lized by a quantum Zeno-like effect. Theoretically, it is 
also possible to modulate periodically the dispersion to 
stabilize the nonautonomous soliton, which is quite read- 
ily realized in the context of optical soliton transmission. 

The above discussions arc entirely based on the 
Painleve integrability condition (2) and are not limited 
to the bright soliton solution. In fact, any allowed solu- 
tion of the standard NLS equation, including the bright 
and dark solitons, the multi-soliton solutions, and even 
the periodic plane wave solutions, can be stabilized by 
such a scheme. This provides a novel approach to study 
in detail the matter-wave properties and related soliton 
dynamics. 

The same scheme can also be applied to the optical 
soliton. The present result has an important implica- 
tion that with the help of the quantum Zeno-like effect 
a dissipationless, ideal optical soliton transmission might 
be possible in reality. Once it is realized, the impact to 
modern information technology is indeed far-reaching. 

The above discussion is based on an intuitive physics 
consideration and the Painleve integrability condition ob- 
tained by a rigorous mathematical derivation. One may 
question its realizability, even be suspicious of it as an ar- 
tifact. In fact, some previous results have an implication 
that our result is physically reasonable. For example, the 
earlier numerical simulations of the nonautonomous NLS 
equation in BEC 27 ' 28 have clearly indicated that the 2D 
bright soliton can be stabilized for some time by tun- 
ing periodically the nonlinearity, even if the instability of 
multidimensional solitons is eventually inevitable 29 . A 
recent experiment in optics 30 also clearly demonstrated 
that the propagation of femtosecond pulses was stabi- 
lized in the layered Kerr media consisting of glass and 
air (the nonlinearity management). Loosely speaking, 
the fact that the dispersion management in nonlinear 
optics can effectively improve the optical soliton trans- 
mission quality is also compatible to the present result. 
It is expectable that an appropriately designed disper- 
sion management can further improve the optical soliton 



transmission quality, and ideally yield a dissipationless 
transmission. 

Method 

Once the Painleve integrability condition equation (2) 
is satisfied, an exact analytical solution of the nonau- 
tonomous NLS equation (1) can be obtained from the 
corresponding canonical solution of the standard NLS 
equation 

i-^Q(X, T)+e-^Q(X, T) + 6 \Q(X, T)| 2 Q(X, T) = 0. 

(4) 

Here e and 6 are constants. When Se > (< 0), the stan- 
dard NLS equation (4) has bright (dark) soliton solutions. 
This is materialized by a general transformation 24 

u(x,t) = Q(X(x,t),T(t))e la ^ t)+c{t \ (5) 

where X(x,t),T(t),a{x,t) and c(t) are real functions 
to be determined by the requirement that u{x, t) and 
Q{X 7 T) satisfy equations (1) and (4), respectively. In- 
serting equation (5) into equation (1) and comparing with 
equation (4), we obtain a set of ordinary differential equa- 
tions for these transformation parameters, whose solu- 
tions read 

°<->^(lKI)-f)* 2 

+Cl 9 j^e^ x - Cle J dt ' 9 j^e^') +C2j{6) 
X(x,t) = ^e^x-2e Cl J dt> f^^'l (7) 
T(t) = Jdt'^e^ + C,, (8) 
c(t) = i In ^|+r(i), (9) 

where T(i) = J Q -f(t')dt' is accumulated dissipation 
and/or gain effect. C\ , Ci , and C3 are arbitrary constants 
and are set to be zero for simplicity. Through this trans- 
formation, all exact solutions, including the canonical 
solitons, of the standard NLS equation can be recast into 
the corresponding solutions of the nonautonomous NLS 
equation (1). This result yields a mapping between the 
canonical and the nonautonomous solitons in a straight- 
forward way. 
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